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ABSTRACT 

Ha emission V/R variations caused by a discontinous mass transfer in interacting binaries 
with a rapidly rotating accreting star are modelled quantitatively for the first time. The program 
ZEUS-MP was used for a non-linear 3-D hydrodynamical modeling of a development of a blob of 
gaseous material injected into an orbit around a star. It resulted in the formation of an elongated 
disk with a slow prograde revolution. The LTE radiative transfer program SHELLSPEC was used 
to calculate the Ha profiles originating in the disk for several phases of its revolution. The profiles 
have the form of a double emission and exhibit V/R and radial velocity variations. However, these 
variations should be a temporal phenomenon since imposing a viscosity in given model would 
lead to a circularization of the disk and fading-out of given variations. 

Subject headings: stars: emission-line, Be - binaries: close - circumstellar matter 



1. Introduction 

The existing modeling of the V/R variations 
was based almost exclusively on the model of 
a slowly revolving elongated envelope around 
a single star, identified in physic al term s with 
one-armed oscillations (see e.g. Struve 1931 
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Johnson! Il958t iMcLaughlinl Il96li lOkazaki 
Fift & Harmaned l2006h in the disk. The 
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lope was assumed to originate from an equatoreal 
mass outflow due to a critical rotation of the star. 
In this paper, we present the first preliminary 
investigation of an alternative idea that an elon- 
gated envelope around a star originates from a 
discontinuous and short-time mass transfer from 
a companion in a binary system. Given transfer 
may occur in eccentric binaries during a periastron 
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passage. Moreover, such an inflow of a material 
could also be caused by a density enhancement 
in a stellar wind from a secondary in a form of 
coronal mass ejections from a solar-like, chromo- 
spherically active secondary. 

In Sect. [21 we present a hydrodynamical mod- 
eling of a discontinuous mass transfer. Radiative 
modeling of an Ha profile orignating in a formed 
disk around an accreting star and a measurement 
of its V / R and radial velocity (RV hereafter) vari- 
ations is presented in Sect. |3l 

2. Hydrodynamical modeling of a discon- 
tinuous mass transfer in a binary 

We decided to carry out the first modeling 
of V/R changes imposed by a discontinous mass 
transfer presented as a blob of a gaseous mate- 
rial put into an orbit around a rapidly rotating 
star having a quadrupole term in its gravitational 
potential. As detailed below, we are using sev- 
eral simplifications which can be challenged but 
we do not think that they seriously affect the re- 
sults on the qualitative basis. Probably the main 
simplification is the assumption of an inviscid gas. 
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A viscosity of an orbiting gas, which is not in- 
cluded in this first model, is expected to destroy 
any asymetry in the disk and all observed effects 
should be only temporal. Another simplification 
which could be criticized is the assumption of an 
optically thin environment. To some defense, we 
would like to mention that in the early st ages of 
the at tempts to model the V/R variations. iHuang 
(1973) used a model based on t he assumpt ion of 
the optically thin envelope while iKfii ( 19761) did a 
similar study assuming optically thick envelopes. 
While the Kffz's line profiles look more realistic, 
both authors obtained a reasonable description of 
the V/R changes since the velocity field and the 
asymmetry of the envelope were decisive ones. On 
the other hand, we should emphasize that we make 
no assumptions about the disk in our attempt and 
gradually build it via a non-linear hydronamical 
modeling of an evolution of a discountinuous mass 
inflow. 

To simulate a hydrodynamical development of a 
blob of mass originated from a discontinous mass 
transfer from a secondary, we used the program 
ZEUS-MP, a multipr ocessor clone of the origina l 
program ZEUS-3D (|Vernaleo fc Reynolds! [20061) . 
A computational domain was a hollow cylinder 
which had the following dimensions in cylindrical 
coordinates {z,r, <p): 

-0.LR* <z< O.LR*, 
R* <r< 70R*, 
< ip < 2tt. 

To its centre, we placed a rapidly rotating star 
with the mass M* = 11 M and the radius i?* = 
5.5 R Q with a gravitational potential <fr(z,r,<p) 
given by 



-GM* 
vV 2 + z 
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where fc 2 denotes an apsidal motion constant and 
/ is a ratio of a surface rotation to a critical Kep- 
lerian rotation: 



/ 



For our simulation, we chose these parameters as 
&2 = 0.03 and / = 0.95. A circumstellar mat- 
ter consists of inviscid, non-selfgravitating ideal 
gas with a constant temperature — 10000 K. 



The dynamical evolution of the gas is described by 
standard hydrodynamical equations for an isother- 
mal case 



= o, 



p— = -Vp-pV$, 

where p,p,v and $ are the gas density, pressure 
and velocity and the gravitational potential. The 
denotes the Lagrangian derivative. Through- 
out the simulation, we used dimensionless scaled 
units (denoted with tildes) for length I, mass m 
and time t. These are related to physical units as 
follows: 




where G denotes the universal gravitational con- 
stant. Note that with this choice we have G = 1 
and the computational domain has the dimensions 
-0.1 < z < 0.1, 1 < f < 70, < <p < 2ir. 

As initial conditions for density and velocity 
components we impose 

p = 10~ 4 , v z =v r = 0, v v = f~ 3/2 

where v v corresponds a keplerian velocity. As 
boundary conditions, we imposed periodic bound- 
ary conditions in tp and z direction. For an in- 
ner boundary condition at the surface of the star 
(f = 1) we took a slip boundary condition v r = 
v z = OjW^ = / with / defined above. Regarding 
an outer boundary condition, we used a pressure 
free outflow boundary condition. 

To simulate the idea of a short discontinuous 
mass inflow to an accreting star, we injected a 
compact (fiaj £ [40,42], <£ inj G [0,0.2], 2 inj G 
[-0.1,0.1]) blob of mass {p = 1) at 80% of the 
Keplerian velocity (v z = v r = 0, v v = 0.8f ~ 3 / 2 ) 
into an initial "interstellar vacuum" and followed 
its evolution in time. 

Our simulation led to a formation of an elon- 
gated accretion disk around the central star. A 
density structure in an equatorial plane of the disk 
and its time evolution is shown on the left side 
of Fig.[T] where are several selected frames from a 
hydrodynamical animation presented in the elec- 
tronic edition. Material in the disk is orbiting 
counter-clockwise and a direction to an observer 
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Fig. 1. — A density structure in an equatorial plane of the disk (left) and the Ha line profiles resulting from 
it (right) for various phases of the prograde rotation of the elongated density enhancement. Gas particles in 
the disk are orbiting counter-clockwise and a direction towards an observer is to the right along a horizontal 
line. See the electronic edition of the Journal for a color version of this figure. The left side figures are also 
available as an mpg animation in the electronic edition. 
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in all panels is to the right along a horizontal line. 
It is seen that a denser region is formed near an 
apocentrum of the elliptical disk due to a crowd- 
ing of the orbiting material there. The elongated 
disk is undergoing a slow prograde rotation with 
a period of about 16 years. 

The prograde rotation of the disk is caused by 
the quadrupole term of the gravitational poten- 
tial. This result is in a qualitative agreement with 
several papers which studied an influence of the 
quadrupole term of the gravitational potencial of 
fast rotating stars o n a precession rate of disks 
around them (see e.g. Savoniie fe Heemskerkll993t 
Okazaki 1996 : Fift fe Harmanec 20061 ). However, 
the linear, one armed oscilation, numerical simula- 
tions presented in these papers use diferent values 
of input parameter and therefore it is not possible 
to do a quantitative comparison of our precesion 
pe riod with their results. M oreover, it was shown 
by iFift fe Harmanec ( 2006 ) that a precesion pe- 
riod is dependent on many free parameters (one 
of them is a term /C2/ which describes a rate of 
stellar distortion due to its rotation) and the de- 
pendence o n some parameters is q uite strong (see 
Figs. 3-9 in lFift fe Harmanedl2006l) . Therefore, by 
a diferent choise of values of given input parame- 
ters, it is possible to obtain a precesion period in 
a range from ^1 year to several tens of years. 

3. Radiative modeling of the Ha profile 

To model the emission profiles of the Ha line 
originating in the disk, we used the program 



SHELLSPEC (|Budai fc Richards! 120041 ) . It solves 
a radiative transfer along the line of sight in an 
optically thin environment assuming LTE. Only a 
scattered light from the central star is taken into 
account, however. The star itself was treated as 
a black body with a temperature T s — 22000 K . 
We chose a view of an observer located in the 
equatorial plane, i.e. an inclination 90°. Result- 
ing profiles of the Ha line for various phases of 
the prograde rotation of the disk are displayed on 
the right side of Fig.Q] The modeling leads to 
double-peaked emission-line profiles. The denser 
region in the disk emits more radiation and its 
revolution around the central star combined with 
the rotation of the disk particles results in the 
V/R changes. 

We imported a representative selection of the 
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Fig. 2. — A time variation of the V/R ratios and 
the RVs for one cycle of the prograde revolution 
of the disk. 



model Ha line pr ofiles into SPEFO (jHorn et al 
1996I : Skoda 1996I ) and measured their V/ R ratios 
and the RVs on wings of the emission. A time 
variation of the V/R ratios and the RVs for one 
cycle of the prograde revolution of the disk is dis- 
played in Fig. [5J It is seen that both quantities 
have a sinusoidal development and that they are 
in phase. 

4. Conclusion 

We tried to test the idea that a temporal Ha 
V/R variations can be caused by a discontinuous 
mass transfer from a companion star in a binary. 
As an initial simple represantation of given dis- 
countinuous mass transfer, we injected a blob of 
mass into an orbit around a rapidly rotating star 
and followed its evolution via 3-D hydrodynamical 
modeling which led to a formation of an elliptical 
disk with a slow prograde revolution. We expect 
that the similar result would be obtained when 
injecting a blob of "new" mass into an already 
developed circular disk. A LTE radiative model- 
ing of the Ha line led to the double-peaked emis- 
sion profile which V/R ratios and RVs undergo a 
cyclic variations in phase. However, these varia- 
tions should be only temporal since a viscosity of 
an orbiting gas should circularize the disk if there 
is an absence of a m ass transfer from a s econdary 
as was also shown bv lBisikalo et al. l <|200lh . There- 
fore also the V/R variations should gradually fade- 
out. 
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Nevertheless, we are aware that this first at- 
tempt to test given idea was considerably sim- 
plified. Beside the already mentioned assump- 
tion of the inviscid isotermal gas, the simplifica- 
tion also involves the way how the discontinu- 
ous mass transfer from the companion was rep- 
resented. Moreover, it is expected that an inclu- 
sion of the attractive force of the orbiting compan- 
ion (also absent in the present model) could speed 
up the disk revolution, resulting in a shorter cycle 
length. Therefore, we intend to improve our mod- 
eling taking into account things mentioned above 
in some future study. 
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